Predicted band structures of III-V semiconductors in wurtzite phase 
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While non-nitride III-V semiconductors typically have a zincblende structure, they may also form 
wurtzite crystals under pressure or when grown as nanowhiskers. This makes electronic structure 
calculation difficult since the band structures of wurtzite III-V semiconductors are poorly character- 
ized. We have calculated the electronic band structure for nine III-V semiconductors in the wurtzite 
phase using transferable empirical pseudopotentials including spin-orbit coupling. We find that all 
the materials have direct gaps. Our results differ significantly from earlier ab initio calculations, 
and where experimental results are available (InP, InAs and GaAs) our calculated band gaps are in 
good agreement. We tabulate energies, effective masses, and linear and cubic Dresselhaus zero-field 
spin-splitting coefficients for the zone-center states. The large zero- field spin-splitting coefficients we 
find may lead to new functionalities for designing devices that manipulate spin degrees of freedom. 

PACS numbers: 71.15.Dx, 71.20.-b, 71.00.00 



I. INTRODUCTION 

Semiconductor nanowhiskers (NW) have attracted a 
tremendous amount of interest in recent year J^^^^^EHEIZlll^ 
Much of it is due to thei r poten tial application in areas 
such as photovoltaic cells^^^^^, nano-electromech anical 
resonator arrayJ^^, m icrowav e and THz detect ioiPlIS, 
single photon detectioi J-^^ * -^^ * -^^ *, field- effectP^^^ and single- 
electron transistors^^J, and variou s other electronic and 
optoelectronic devicej^ ^^ ' ^^ ' ^^ ' ^^l NWs are also inter- 
esting because ID systems can be made using material 
combinations for which large lattice mismatches prohibit 
quantum well structures, allowing greater freedom in ma- 
terial combinations for device engineering. 

In contrast to bulk non-nitride III-V materials which 
are usually zincblende (ZB), NWs predominantly crystal- 
ize in the wurtzite (WZ) phase^^^^^^. Several different 
theoretical explanations for this behavior have been pro- 
posed. Recent calculations suggest that the WZ phase 
is energetically favorable for small NW radii^^^^, al- 
though this does not account for all NW radii for which 
the WZ phase is experimentally observed. Calculations 
based on an empirical nucleation model indicate that WZ 
formation is favored for certain ranges of the interface 
energie^^. Ab initio calculations indicate that the WZ 
phase is favored due to the accumulation of electrons at 
the interstitial site containing the Au catalyst^^, while 
other calculations show that the polytype is determined 
by growth kinetics^^. It should be noted that these differ- 
ent mechanisms are not necessarily mutually exclusive. 

Theoretical understanding of the electronic and optical 
properties of semiconductor nanostructures is based on 
a knowledge of the electronic properties of bulk materi- 
als. However little is known about the electronic band 
structure of III-V semiconductors in WZ phase since 
most do not naturally occur as bulk crystals. Moreover, 
the NWs often contain sections of ZB material, form- 
ing heterostructures out of the differing band structures 
of the two polytypes ^ ' ^^'^^ ' ^^ . Band structures of WZ 



III-V semiconductors have been calculated using density 
functional theory (DFT) in the local density approxima- 
tion (LDA^P^f^ Since the LDA underestimates band 
gaps the WZ band structure can not be directly deter- 
mined, and instead calculations of WZ and ZB are typi- 
cally compared to obtain the differences between the two 
polytypes. The band structures of GaAs and InAs in WZ 
phase have also been calculated using the GW methocPl, 
giving somewhat different results than those from the 
LDA. In addition to the inherent errors in ab initio band 
gaps, all of the above calculations neglected the spin- 
orbit coupling, which is known to significantly alter the 
valence band structure of semiconductors. 

In this paper we present calculations of the bulk elec- 
tronic band structures of the nine non-nitride III-V semi- 
conductors in WZ phase using empirical pseudopoten- 
tials including spin-orbit coupling. These calculations 
are based on transferable model pseudopotentials assum- 
ing ideal WZ structure. The spherically symmetric ionic 
model potentials are first obtained by fitting the calcu- 
lated bulk ZB energies to experimental energies at high 
symmetry points. The band structure of the WZ poly- 
type is then obtained by transferring the model pseu- 
dopotentials to the WZ pseudopotential Hamiltonian us- 
ing the appropriate crystal structure factors. 

This method has been proven to be very successful 
in obtaining the bulk band structures of semiconductor 

polytype&40 41 42 43 44 45 46 47 48^ rj^^^ ^^^^^ ^^^^^^ 

pseudopotentials are specific to each material and are 
only transfered between polytypes. Therefore the model 
potentials should be transferable between ZB and WZ 
polytypes due to the similarities in their crystal struc- 
tures. In both structures all of the nearest neighbors and 
nine out of the twelve second nearest neighbors are at 
identical crystallographic locations^^ while the the sec- 
ond nearest neighbors are equidistant. 

This paper is organized as follows. In section [n| we out- 
line the similarities and differences between ZB and WZ 
crystal structures as well as the direct correspondence be- 
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tween high symmetry /c-points in the two poly types. In 



section |ni] we describe the transferable pseudopotential 
method. In section IV we present the calculated band 
structures, their respective density of states (DOS), and 
effective masses. Finally, we summarize the results in 
section [Vll 



II. WURTZITE VS ZINCBLENDE 

A. Crystal Structure 

The ZB crystal is formed by two interpenetrating face 
centered cubic (FCC) Bravais lattices (each of a different 
atomic species) , whereas the WZ structure is constructed 
from two interpenetrating hexagonal close-packed (HCP) 
lattices. The differences between the two structures 
are best understood by viewing along the [111] direc- 
tion (Fig. [l^,b), along which both look like stacked 
hexagonal layers. The atoms are identical within each 
layer, and the layers alternate between the anion and 
the cation. For the ideal WZ crystal the lattice con- 
stant is given by awz = clzb/ and the lattice constant 
along the c-axis (axis perpendicular to the hexagon) is re- 
lated to the in-plane lattice constant by c = (8/3) awz- 
Since the WZ crystal is tetrahedral, the nearest neigh- 
bors are the same in the two polytypes. In addition, 
we see from Fig. [l] that nine of the twelve nearest 
neighbors in WZ are the same as in ZB. These struc- 
tural similarities suggest that the local electronic en- 
vironment will be the same in the two crystals, and 
therefore the crystal potentials will be nearly identical 
in WZ and ZB. In WZ the type-1 atoms are located at 
(0,0,0) and |ai + + ^as, while the type-2 atoms 
are located at usi^ and |ai + |a2 + {\ + 'u)a3, where 
the primitive lattice vectors are ai = (1, V^, 0)avt^^/2, 
a2 = (1, — \/3, 0)a\yz/2 and as = (0,0,c). The param- 
eter u = 3/8 for the ideal WZ crystal structure, which 
we will assume throughout this paper. This assumption 
is supported by the fact that experiments show that WZ 
GaAs is very close to ideal^^. 



B. Band Structure 

Due to the similarities of the two crystals, many of the 
high symmetry points in the Brillouin zones of ZB and 
WZ are related to each other, and an understanding of 
their correspondences is useful for understanding trends 
in the band structures of the two polytypes. Fig. |2]shows 
the relationships among the zone center states in WZ and 
the corresponding L and F points in ZB, both with and 
without spin-orbit coupling. One of the most important 
features is that in the empty lattice approximation the 
L-point in ZB is zone- folded to the F-point in WZ. As 
a result, in the absence of spin-orbit coupling, the Fi, 
Li and L3 states in ZB correspond to Fi, F3 and F5 
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FIG. 1: (a) Staggered configurations of atom-1, as 
viewed along [111] for ZB (b) Eclipsed configuration of 
atom-1 as viewed along [111] for WZ. Note that 9 out of 

the 12 second nearest neighbors are in the same 
position. The other 3 are rotated by 7r/3. (c) Brillouin 
zone for ZB. (d) Brillouin zone for WZ. F is at the 
center of the Brillouin zones. 



respectively in WZ (Fig. [2|. Because of this folding over 
of the L- valley, indirect gap ZB materials with an L valley 
conduction band minimum would be expected to have a 
direct gap in the ZB phase unless the energy of the state 
was significantly shifted by the crystal potential. 

The states at the top of the valence band in WZ also 
have some important differences with their ZB counter- 
parts. In the absence of spin-orbit coupling, the hexag- 
onal crystal field of WZ splits the p-like Fis state of ZB 
into a four-fold degenerate Tq and a doubly degenerate 
Fi. In terms of the p-orbitals, these states are Pz Ti 
are Px^Py Fg. With the inclusion of spin-orbit cou- 
pling, splits into the Fg^ heavy-hole and the T^v light- 
hole. Therefore, all zone center states in WZ belong to 
either F7, Fg, or Fg. There are similar correspondences 
between the high symmetry directions of the two crys- 
tals. The symmetry line A(F L) in ZB corresponds to 
the A(F ^ A) line in WZ^. 
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FIG. 3: (a) Zero field Dresselhaus spin splitting in WZ 
along kx for a F7 conduction band, (b) The zero field 
spin splitting is not seen along the kz direction. 

(F A), the crystallographic point group is Cq^^^. As 
all the irreducible representations compatible with spin 
(F7^8,9) in this group are doubly-degenerate^^, there is no 
spin-splitting along kz- 

In WZ, the spin splitting effects near zone cente r can be 
described with an effective one band HamiltoniaiPI^^IIIl 



C. Spin Splitting 

For certain individual states in crystals lacking inver- 
sion symmetry, spin-orbit coupling causes a splitting of 
spin- up and spin-down states which leads to E^{—\^) ^ 
^l(k) ( see Fig. [3|. The states at k = are still two- 
fold degenerate, resulting in a nonvanishing Vfc£^(k) at 
the origin for certain crystallographic directions. In the 
case of WZ, the states remain spin-degenerate for k along 
the c-axis because at any point along the kz direction 



with SO ^^^^^gfipflfl^^ ^jth SO 
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FIG. 2: Schematic showing the correspondence between 
energy levels at the L and F points in ZB and the F 
point in WZ with and without spin-orbit coupling^''. 
The dashed lines show the correspondence between the 
states. Additional degenerate levels are shown in grey. 



i/(k) oc ^^ + ^ + ^a4a.fc,+cT,fcj(2"+i) (1) 
2m^ 2m\\ ^-^ 

where are constants. This effective Hamiltonian is 
invariant with respect to Cq^ for all odd powers of k. 
The coefficients of the linear and cubic Dresselhaus spin- 
splitting terms can be obtained by expanding Eq. [l]up 
to n = 1. Near the crossing, the difference in energies 
between the spin-split bands in Fig. (sj^a) has the fornP^ 

iE,,t(fc)-^u(fc)i = 2ciA: + dfc' (2) 

where ("i and Q are the linear and cubic Dresselhaus 
spin-splitting coefficients respectively for band i. These 
coefficient characterize the spin spin splitting, and may 
be determined from the computed band structure by 
fitting a function of the above form. While the exis- 
tence of a spin splitting in a particular direction is de- 
termined by symmetry consideration^^, the magnitude 
can vary considerably. Although both ZB and WZ crys- 
tals lack inversion symmetry, spin splitting effects are 
much mo re prominent in WZ due to its lower crystal 
symmetr^™^ ^^l^^ l ^^ l ^" l ^^ l ^'^ l ^^ l ^^'^^ l ^^' , which can be of 
particular use in spintronic devices. 

III. METHOD 

A. Pseudopotentials 

Our band structures are computed using the empiri- 
cal pseudopotential method of Cohen and Chelikowsk}!^ 
with a model potential applicable to both ZB and WZ 
structures. Pseudopotentials exploit the fact that the 
electronic wave function may be separated into the sum 
of a rapidly oscillating part near the atomic cores and a 
slowly varying piece. The pseudopotential approach re- 
lies on the assumption that the core electrons are frozen 
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and that the valence electrons move in a weak single- 
electron potential making the true atomic wave function 
orthogonal to the core states. The pseudo wave equation 
is then 

^2 



P 

2m 



Vr 



pp 



El 



(3) 



where is the smoothly varying pseudo wave function, 
and Vpp is the pseudopotential which includes a repulsive 
core to partially cancel the deep potential near the atomic 
cor^^. By taking Vpp to be a local pseudopotential which 
is a function only of position, it may be expanded in terms 
of reciprocal lattice vectors, G, as 



V(r) = Y.Vj^{G)S^{G)e' 



Gr 



G,a 



1 = 1 



(4) 



(5) 



where a labels the atom type, V^^ (G) is the form factor, 
Sa{G) is the structure factor, is the number of atoms 
per unit cell of type a, and r^^i is the position of atom 
number i of type a. 

For binary compounds it is convenient to separate the 
pseudopotential into symmetric (S) and antisymmetric 
(A) parts as 

{G'\Vpp\G) = Vs{G' -G)Ss{G' -G) 

+ iVA{G^ -G)Sa{G' -G). (6) 

The symmetric and antisymmetric structure factors are 
given by 

Ss{G) = ^y^expHG-r,) (7) 



N ^ 
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(8) 



where N is the number of atoms per unit cell and Pj = 
+1 for one type of atom and —1 for the other type. The 
symmetric and antisymmetric form factors, Va(G) and 
Vs(G\ are obtained from the sum and difference of the 
spherically symmetric anion and cation potentials. 

There are several approaches to calculating Vp^^^. 
In the empirical pseudopotential approach used here, 
VsiG) and Va(G) are adjusted to fit the calculated en- 
ergy spectrum to experimentally determined energies at 
band extrema. Model potentials that yield an accurate 
band structure of a known polytype should reliably pre- 
dict the band structure for the unknown polytype if the 
two crystal structures are similar. To compute the WZ 
band structure using pseudopotentials obtained from ZB 
requires Ks'(G) and Va{G) to be continuous functions 
that can be evaluated at any value of G. A w ide variety 
model potentials have been used ^^ ' ^^ ' ^^ * ^^ * '^-'^* and we use 
potentials of the form 

ys(G) = (9) 



Va{G) 



exp(x3G^ + X4) + 1 
{x[G'^ + 4) exp(4G^ + X4) (10) 



where G = |G|, and the parameters Xj and Xj are ob- 
tained for each material by fitting to the ZB band struc- 
ture. 



B. Spin-Orbit Interactions 



We include the spin-orbit coupling, given by 



(W(r) xp)-a 
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(11) 



cannot be used 



using the method of Weisz^^. Eq. 
directly since spin-orbit coupling involve the core states 
which are omitted from the pseudo wave equation as a 
result of orthogonalization. Instead, by returning to the 
original Schrodinger equatio n we m ay expand in terms of 
the core states |(/)c) to obtairP^^ 



b'\Vso\<l}) - Y.^4>'\4>c'){4>AHso\4>c){cl>c 



(12) 



In this way the matrix elements of Hso for core states may 
be parameterized and then fit to experiment in the same 
way as the pseudopotential form factors. By expanding 
1^) in Bloch functions, and expressing the core states in 
terms of its constituent atomic radial wave functions and 
respective spherical harmonics, Eq. 12 can be recast a^ 



{'K',s'\Vso\K,s) = (K' X K) • (^Vl^) 

^Xi P/(cos^k'.k)^(K'-K) 

I 

(13) 

where P/ is the derivative of the Legendre polynomial, 
crs are the Pauli matrices, K = G + k, and is the angle 
between K and K^ The coefficient is given in terms 
of the core wave functions by 



Xi=fi(3ni{K^)/3ni{K) 



Pnl{K)=C 



/»oo 

Jo 



(14) 



l)jni{Kr)Rni{rydr (15) 



where C is a normalization constant such that /3ni (K) / K 
approaches unity in the limit K goes to zero and n is 
the principal quantum number for the core state being 
considered. For III-V semiconductors, it is not required 
to expand Eq. [13] beyond / = 2 since they do not have 



core shells filled beyond d-orbitals. Expanding Eq. 13 up 
to / = 2, the matrix elements for the spin-orbit coupling 
in a binary compound are 



{'K',s'\Vso\K,s) = X K) . {s'\ct\s) 

(A,^ + AfK'-K) Ss{G'-G) 

(a^ + A^K^-k)5^(G'-G)" 



(16) 
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(AW + Ap))/2 

(AW-Ap))/2 

W/?iV(K.)/3iV(K,) 
7, w/?i^)(K,)/3i^)(K,) 



(17) 

(18) 

(19) 
(20) 



where the superscript (1),(2) specifies which atom, the 
coefficient jui is an empirically adjusted parameter and 
7/ is the ratio of the anion to cation spin-orbit splitting 
energies for a given core statd^. The overlap integral, 
Pnii is constructed from the atomic core wave functions 
The radial part of the core wave function. 
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using Eq. 

Rni, is an approximate Hartree-Fock solution taken from 
Herman- Skillman tables^^. For Ga, In, As and Sb, terms 
up to / = 2 in Eq. [16] are included while Al and P only 
go up to / = 1 since they do not have valence d-shells. 



With the inclusion of Eq. [16] the total pseudopotential 
Hamiltonian becomes 



H 



2m 



Vpp ■ 



(21) 



C. Fitting 

The pseudopotential parameters x^, and /i^ in 
Eqs.[9]and[l0 were determined by fitting the band struc- 
ture obtained from the Hamiltonian of Eq. 21 to exper- 
imental energies of the band extrema of ZB materials. 
The Hamiltonian was evaluated in a plane wave basis 
with a cutoff of |G| < 327r/a, and for each value of k 
the Hamiltonian was diagonalized to give energies to be 
fit to the experimental target values. The fitting was 
accomplished by minimizing the error function 



[Ei (calculated) — Ei{target)Y 
Ef (target) 



(22) 



where the sum over i ranges over the targeted en- 
ergies Ei, and Wi are weighting factors adjusted to 
speed convergence, 
ting parameters, (all 



z?r TT^r Tj^x Tj^L TT^r TT^r T^r 

^6c' ^7v^ ^6c^ ^6c' ^v6^ ^c7^ ^c8' 



Seven energies were used as fit- 
with respect to E^^ = 0): 
The targeted values of 



Eq^, Eq^, Eq^ and E^y were taken from Ref. 76] while the 
higher transition energies, were taken from Ref. [77l In 
addition constraints were imposed to ensure the correct 
band ordering of valence states by forcing the third and 
fourth (spin-degenerate) valence band states to have Fy 
symmetry. F was minimized with respect to x^, and 
jiii using Powell's method, with the local pseudopotential 
form factors from Ref. 78 used as an initial starting point. 
Slightly different initial values were used as a check that 
the solution did not converge to a spurious local mini- 
mum since Powell's method finds the local minimum. In 
those cases where the value of F was large, indicating a 
poor fit, the weights Wj were adjusted to climb out of 
the local minimum in which the algorithm was trapped, 
and the minimization algorithm was continued. 



D. Transferable Pseudopotentials 

Once the form factors have been determined for the ZB 
poly type, they may be transfered to the WZ structure by 
centering the spherically symmetric atomic pseudopoten- 
tials on the positions of the ions in the WZ form. The 
transferability depends on the similarities of the crystal 
structures, and for sufficiently dissimilar polytypes one 
would expect the method to fail. Fortunately, as dis- 



cussed in section II A ZB and WZ are very similar, as the 
WZ crystal structure can be thought of as a variation of 
ZB with the same local structure, but a slightly different 
long range structure. It is important to note that all the 
parameters were fit independently for each material, and 
the pseudopotentials were transferred between polytypes 
with the same binary composition. Since the WZ prim- 
itive cell has four atoms (rather than the two of ZB), 
its structure factor contain more terms. Substituting the 
atomic positions in section [ITA] into Eqs.[7]and[8]the WZ 
structure factors are 



Ss = \{l + e- 
SA = \{l-e 



-iGsuc 
-iGsuc 



)(^l^e-'G,a/VS+Gsc/2^^ (23) 
) (^l + g-^G'2a/V3+G3c/2^ ^24) 



where Gj {j = 1, 2, 3) are the components of the recipro- 
cal lattice vector G. 



IV. RESULTS 

A. Calculated III-V Zinc Blende band structures 

The pseudopotential parameters determined by fitting 
to the zincblende band structures are given in Table [T] 
including the spin-orbit parameters. The accuracy of the 
results may be gauged by Table [Til which gives the ratio 
of each band energy to the experimental value to which it 
was fit^. We see that the results agree with experiment 
to within 1% for all but Eq^. Even for Eq^, the deviation 
from the experimental value is greater than 10% only for 
InSb. We have also made a similar comparison of the ef- 
fective masses and Luttinger parameters'^. These values 
differ from the energy ratios in that the masses were not 
fit to experimental data. The masses were determined 
by doing a quadratic fit to the band extremum, with the 
Luttinger parameters determined from 



[001] 



7i T 272 



'^hhiih) J 
* \ [111] 

= 71 ± 273 

^hh{ih) J 



(25) 



It should be noted that several of the masses listed in Ref. 
[76l are either obtained theoretically or have large experi- 
mental uncertainties. For example, only theoretically cal- 
culated effective masses are available for AlP and for the 
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valence band of AlSb. To the best of our knowledge, ex- 
perimental results are not available for these compounds. 
In the case of GaP the conduction band effective mass 
is extrapolated from its ternary alloy. For InAs, there 
is great experimental uncertainty about the heavy and 
light hole masses^^. Even though we have omitted non- 
local corrections not associated with spin-orbit coupling, 
our effective masses are in very good agreement for com- 
pounds containing lighter elements. 



B. Predicted III-V Wurtzite band structures 

The calculated band structure and the corresponding 
density of states(DOS) for each of the nine III-V semicon- 
ductors in WZ phase are shown in Fig. [4][l2] The elec- 
tronic band structures are calculated in the irreducible 
wedge of the Brillouin zone (Fig. [ijd)). It should be 
noted that the band structure of WZ is more compli- 
cated than that of ZB due to it lower crystal symmetry 
and has roughly twice as many bands over a given energy 
range. The irreducible representations of the zone center 
states was determined by transforming the pseudo wave 
functions under the symmetry operations of the respec- 
tive crystallographic point group. 

Tables |III|[V| list the calculated zone center energies 
of the band extrema, effective masses for k parallel 
and perpendicular to the c axis, and the linear and 
cubic Dresselhaus spin splitting coefficients and ^2- 
We calculated the Dresselhaus coefficients by fitting a 
function of the above form in Eq. |2] to the calculated 
band structures. These parameters may be used in 
constructing k-p WZ Hamiltonians for nanostructure 
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FIG. 4: (a) Calculated low temperature band structure 
for AlP in WZ phase, (b) Calculated DOS. 
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FIG. 5: (a) Calculated low temperature band structure 
for AlAs in WZ phase, (b) Calculated DOS. 
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FIG. 6: (a) Calculated low temperature band structure 
for AlSb in WZ phase, (b) Calculated DOS. 



calculationj^ ^ ' ^" ' ^^ ' ^^ ' ^^ ' ^^ ' ^^ !. We have also complied a 
summary (Table |VI[ ) listing the energy differences most 
important for nanostructures. 

Foremost is the band gap and the irreducible repre- 
sentation of the conduction band minimum. In our cal- 
culations it is seen that all of the materials containing 
Al or Ga have Fg conduction band minima in the WZ 
phase, whereas all of the materials containing In have F7 
conduction band minima. We also give the spin-orbit en- 
ergy, A50, the crystal field splitting, Acr, and the offset 
between the valence band edge of each polytypes. A50 
and Acr are extracted using the quasi cubic approxima- 
tion which assumes WZ to be equivalent to [111] -strained 
2£79 86^ with Ago and Acr related to the Vjy hole ener- 
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Material 


Xi 


X2 


X3 


X4 






4 


X4 


Ml 


/i2 


AlP 


0.083 


-0.579 


0.031 


-2.586 


-0.11 


0.9 


-0.061 


-1.178 


0.012 





AlAs 


0.062 


-0.459 


0.027 


-2.629 


-0.041 


1.003 


-0.056 


-1.693 5.28 X 10"^ 


5 X 10"^ 


AlSb 


0.06 


-0.412 


0.032 


-2.548 


-0.09 


0.17 


-0.051 


-2.098 


7.26 X 10"^ 


2.96 X 10"^ 


GaP 


0.085 


-0.457 


0.04 


-2.566 


-0.351 


5.165 


-0.205 


0.339 


0.385 


-2.3 X 10"^ 


GaAs 


0.058 


-0.467 


0.023 


-2.583 


-0.063 


1.091 


-0.074 


-1.298 


0.052 


8.3 X 10"^ 


GaSb 


0.042 


-0.343 


0.022 


-2.584 


-0.009 


0.618 


-0.043 


-2.233 


0.056 


2.78 X 10"^ 


InP 


0.049 


-0.385 


0.027 


-2.602 





0.847 


-0.059 


-1.654 


0.243 


-1.09 X 10"^ 


InAs 


0.036 


-0.298 


0.033 


-2.615 


-0.011 


1.359 


-0.121 


-1.124 


0.082 


2.603 X 10"^ 


InSb 


0.022 


-0.174 


0.023 


-2.42 


-0.012 


1.158 


-0.082 


-1.363 


0.085 


5.7 X 10"^ 



TABLE I: Fitting parameters for symmetric and antisymmetric form factors, where the form factors are in units of 
Ry. jii and 112 are the fitting parameters for the spin-orbit coupling. All fitting coefficients are rounded off to the 

third decimal place. 



Material 



AlP 
AlAs 
AlSb 
GaP 
GaAs 
GaSb 
InP 
InAs 
InSb 



IT 



^8c 



1.00 1.00 

1.00 1.00 

1.00 1.00 

1.00 1.00 

1.00 1.00 

1.00 1.00 

1.00 1.00 

1.00 1.00 

1.00 1.00 



1.00 1.00 

1.01 1.00 
1.00 1.00 
1.00 1.00 



1.00 
1.01 
1.00 
1.00 
1.01 



1.00 
1.00 
1.00 
1.01 



1.01 
1.00 
1.00 
1.00 
1.00 
1.00 
1.00 



1.00 1.00 



1.01 
1.00 
1.00 
1.00 
1.00 
1.00 
1.00 
1.00 



1.00 
1.00 
1.00 
1.00 
1.00 
1.00 
1.00 
1.00 



1.03 
0.99 
0.97 
0.97 
0.95 
0.94 
1.02 
0.95 
0.84 



71 72 73 



0.72 
0.99 
1.04 
1.00 
1.16 
1.47 
1.08 
1.36 
1.85 



1.06 0.92 

1.07 0.95 
1.30 0.67 

1.02 0.96 
1.19 0.81 
1.56 0.63 

1.03 0.96 
1.05 0.61 
1.73 0.55 



0.99 0.95 
0.84 0.96 
0.66 0.70 
1.66 1.20 
0.90 0.83 
0.61 0.62 
0.98 1.00 
0.60 0.62 
0.54 0.56 



TABLE II: Ratios of calculated/targeted band transition energies at various high symmetry points for nine III-V 
zincblende semiconductors. All transition energies are referenced to the top of valence band, (Eg^). The spin-orbit 
energies are Ago = E^^ — E^y and A^^ = E^^ — Ej^. In addition, the effective for the conduction band (mc), split-off 
bands(m5o), and heavy and light holes are compared to those from Ref. 76, The heavy and light hole masses are 
expressed in terms of Luttinger parameters 71,72 and 73 (see Eq. 25). Note, that the effective masses were not set 

as targets for fitting the form factors. 
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7: (a) Calculated low temperature band structure 
for GaP in WZ phase, (b) Calculated DOS. 
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FIG. 8: (a) Calculated low temperature band structure 
for GaAs in WZ phase, (b) Calculated DOS. 
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FIG. 9: (a) Calculated low temperature band structure 
for GaSb in WZ phase, (b) Calculated DOS. 



FIG. 11: (a) Calculated low temperature band structure 
for InAs in wurtzite phase, (b) Calculated DOS. 
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FIG. 10: (a) Calculated low temperature band structure 
for InP in wurtzite phase, (b) Calculated DOS. 



FIG. 12: (a) Calculated low temperature band structure 
for InSb in wurtzite phase, (b) Calculated DOS. 



gies by 



^(F^;^)-^(F9.) 



± 



^so'-^cr 



(26) 



where is a parameter relating the lattice constants, 
^/u = a/c. For the ideal WZ structure assumed here, 
and is 3/8 for an ideal WZ structure. 

In ZB phase, AlP, AlAs, AlSb and GaP are indirect 
gap semiconductors with their conduction band minima 
ordered X, L, F. Our calculations show that in all of 
the indirect gap ZB semiconductors become direct gap 
WZ materials with Fg conduction band minima. Previ- 
ous LDA calculations obtained indirect gaps for AlP and 



AlAs in W Z phase, with the conduction band minimum 
at Although this is not the case in our results, 

we do find the M valley conduction band minimum only 
slightly above the F minima (182 meV for AlP and 157 
meV for AlAs). The same LDA calculations also pre- 
dicted AlSb and GaP to have Fg conduction band min- 
ima, in agreement with our results^^. The direct gaps 
of AlP and AlSb in WZ phase are larger than the indi- 
rect (X-valley) gaps of their respective ZB counterparts. 
Whereas the direct (Fg) gaps of AlAs and GaP in WZ 
phase are smaller than the indirect gaps of their ZB coun- 
terparts. 

Although GaAs and GaSb have direct gaps in the ZB 
phase, for both materials the zone- folded L valley in WZ 
is lower, giving a Fg minimum. For WZ, GaAs has a 



band gap 1.5 eV) which is only shghtly smaher than 
that of its ZB polytype, while GaSb has a significantly 
smaller gap (~ 0.53 eV) in WZ phase. In both cases 
this is due to £^lic being close to Eric- This behavior is 
most apparent for GaSb, hence our prediction for GaSb 
is consistent with that of Ref. |38] and Ref. |37l In light 
of recent experimental results, the case of WZ GaAs will 
be discussed in greater detail below. All three Indium 
containing compounds InP, InAs and InSb are direct gap 
semiconductors in WZ phase with Ft conduction band 
minima and have higher band gaps than their respective 
direct gap ZB polytypes. 

While little is known about the WZ polytypes, experi- 
mental data is available for the band gaps of GaAs, InAs 
and InP. Table ( VII ) compares our calculated band gaps 



for these three materials to experiment and ab initio cal- 
culations. We have only tabulated their SX and GW 
results as they make different qualitative predictions as 
to whether the band gap of WZ GaAs is larger or smaller 
than for ZB. Note that the values listed for the ab initio 
methods are not listed directly in the respective refer- 
ences. As none of the directly obtained gaps from ab 
initio calculations are correct, the values are obtained 
from the experimental ZB gap and the calculated per- 
centage change between the two polytypes. As can be 
seen, our results are in very good agreement with experi- 
ments. More importantly, they are in agreement with the 
experimental trends as to whether the WZ band gap is 
larger or smaller than the ZB gap. This is most apparent 
in the case of GaAs, for which ab initio calculations (ex- 
cept the SX methocP^ predict that GaAs should have 
a larger direct band gap in WZ phase than in ZB, in 
disagreement with experiment. All of the low temper- 
ature experimental results show that for GaAs the WZ 
phase has a smaller band gap than ZB. In addition, the 
experimental gaps are obtained from photoluminescence 
measurements on GaAs nanowires, which will be slightly 
larger due to confinement. 

In the case of InP our results are in agreement with 
experiments as well as with the trends from LDA cal- 
culations. Photocurrent spectroscopy measurements on 
InAsi_xPx nanowires have been extrapolated in x to ob- 
tain a band gap for InP of 1.645 eV^, which is higher 
than valules from photoluminescence measurements. The 
fact that the measured values are higher than our pseu- 
dopotential calculations may be du to confinement effects 
in the nanowire experiments. Our calculated gap for InP 
is also in good agreement with the calculations of Zanolli 
et al^\ 



The spin-orbit coupling can alter the ordering of the 
valence band states in different WZ semiconductors. For 
example, in CdS and CdSe the top three valence states 
are in descending order Fg, F7, F7 (referred to as nor- 
mal ordering)^^^ while in ZnO the ordering is F7, Fg, 
F7 (anomalous ordering) which results from a negative 
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E, (eV) 


A,o (eV) 


Aer (eV) 


/\EcB (eV) 


AlP 


2.969 (Pg) 


0.070 


0.409 


-0.1428 


AlAs 


1.971 (Pg) 


0.319 


0.338 


-0.0841 


AlSb 


1.891 (Pg) 


0.683 


0.276 


-0.0909 


GaP 


2.251 (Pg) 


0.082 


0.341 


-0.0801 


GaAs 


1.503 (Pg) 


0.351 


0.244 


-0.0632 


GaSb 


0.509 (Pg) 


0.777 


0.239 


-0.1186 


InP 


1.474 (Pr) 


0.108 


0.303 


-0.0646 


InAs 


0.481 (Pr) 


0.379 


0.195 


0.0405 


InSb 


0.287 (Pt) 


0.787 


0.159 


0.0872 



TABLE VI: Energies of the direct gap III-V WZ 
semiconductors. The symmetry of the conduction band 
minimum is indicated with the band gap. A^o and Acr 
are the spin-orbit sphtting and crystal-field splitting 
energies extracted using Eq. (26). 
/S.Ey B = — Ey"^ , is the energy difference between 
the top of the valance bands for the two poly types. 



spin-orbit energy. In our calculations, all materials ex- 
cept InSb have normal ordering. In the case of InSb the 
ordering of the valence band states is complicated by its 
very large spin-orbit splitting which forces the Pr split- 
off hole bellow the next Pg state (which comes from the 
'folded over' p-\ike L- valley states in ZB). This results 
in the unusual Pg, P7, Pg, P7 ordering of valence band 
states in InSb. 

Tables |III||V| also give the Dresselhaus coefficients for 
each band. As expected^^, all Pg states have zero linear 
Dresselhaus coefficients (Ci) while all Py and Pg states 
have nonzero linear and cubic spin-splitting coefficients. 
Amongst the five WZ phase semiconductors with the 
largest spin-orbit energies (A50), InSb has the largest 
Dresselhaus coefficients, followed by GaSb, InAs, GaAs 
and AlSb. 



cal pseudopotentials with the inclusion of spin-orbit cou- 
pling. The predicted band structures are based on the 
concept of transferable model potentials. Our calcula- 
tions show that in WZ phase, InP, InAs and InSb have 
a direct gap (P7) which is larger than the corresponding 
zincblende material. AlP, AlAs, AlSb, GaP, GaAs and 
GaSb also have direct gaps, but with Pg conduction band 
minima. WZ AlP and AlAs have larger direct gaps than 
the indirect gaps of their ZB polytypes. The opposite 
trend is seen in AlSb and GaP which have smaller direct 
gaps than their indirect gaps in ZB phase. In WZ phase 
GaAs and GaSb have direct gaps which are smaller than 
their ZB counterparts. 

Our calculations are in excellent agreement with with 
experimentally obtained band gaps for GaAs, InAs and 
InP in WZ phase. Significantly, our results agree with 
experiment over whether the WZ band gap is larger or 
smaller than the ZB gap, in contrast to ab initio methods. 
We have extracted the linear and cubic Dresselhaus spin- 
splitting coefficients and find they are generally, though 
not always, larger for materials with larger spin-orbit 
coupling. The relatively large spin-splittings may be 
of use for spin-dependent transport in WZ nanowires. 
More recently, the WZ phase of GaAs has been grown in 
bulk^^. More experimental measurements on such bulk 
WZ phase III-V semiconductors, would lead to a clearer 
understanding of their electronic properties. 



SUMMARY 
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